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ABSTRACT 

Very massive stars are thought to be formed in the early Universe because of a lack 
of cooling process by heavy elements, and might have been responsible for the later 
evolution of the Universe. We had an interest in vibrational stability of their evolution 
and carried out the linear nonadiabatic analysis of radial and nonradial oscillations for 
population III very massive main-sequence stars with 500 — 3000Mq. We found that 
only the radial fundamental mode becomes unstable due to the e-mechanism for these 
stars. The instability appears just after the CNO cycle is activated and the nuclear 
energy generation rate becomes large enough to stop the pre-main-sequence contrac- 
tion, and continues during the early stage of the core hydrogen burning. Besides, we 
roughly estimated amount of mass loss due to the instability to know its significance. 
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1 INTRODUCTION 

The first-generation stars born in the early Universe are 
thought to have no, or few heavy elements. They might play 
an important role for chemical evolution of the early Uni- 
verse. The initial mass function for the first generation stars 
has not been definitively determined yet. However, a lack of 
heavy elements leads to a deficiency of the coolants in the 
star-formation stage and hence to an expecta tion that very 
mass i ve stars might have been formed (e.g. iBromm et al.l 
1 19991 : lAbel et al. I l2002l : lOmukai fc Palld 12003 ). 

iLedouxl |l94ll ) suggested that very massive stars might 
experience vibrational instability in the core hydrogen- 
burning stage. In case of very massive stars, influence of 
radiation pressure makes the effective ratio of specific heat, 
7, close to 4/3. Then pulsation amplitude in the core be- 
comes comparable with that in the envelope. This situation 
is favorable for instability due to the e-mechanism since ex- 
citation by this mechanism in the core can exceed flux dis- 
sipation in the envelope. Many authors have analyzed the 
stability of very massive stars and proposed the critical mass 
above which the instability occurs. With the pure electron 
scattering opacity the radial f undamental mode becomes un - 
stable in case of > 6OM jSchwarzschild fc Harmlll959h 
while with the recent OPAL opacity table ( Rogers fc Iglesiaa 
1 19921) the criti cal mass is 121Mq for the solar metallicity 
(Stothcrs 1992). In the latter case, the fundamental mode is 
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excited mainly by the K-mechanism of the metallic absorp- 
tion rather than the e-mechanism. 



In the zero-metalli city case, how e ver, t he metallic ab- 
sorption never occurs. [Baraffe et al.l (J2001 ) analyzed the 
stability of population III very massive stars and found 
that the fundamental mode becomes unstable due to the 
e-mechanism for stars with 120 — 500Mq. 



Population III very massive stars might have released 
a significant amount of heavy elements by supernova ex- 
plosions and contributed to the chemical evolution of the 
early Universe. The stars with > 300Mq experience the 
core-collapse supernova explosion at the end of their life, 
while those wi th 130 — 300M,^ th e pair-instability super- 
nova explosion (|Ohkubo et al.ll2006l ). In addition, the former 
stars leave intermediate mass black holes, for which several 
candidates have bee n found as ultraluminous X-ray sources 
jFeng fc Soriall201llh Then, the evolution of the very mas- 
sive stars is very important for the above phe nomena, and 
have been investigated b y several groups (e . g. iKl app 1983, 
19841 : iHeger et al'l ,2003 1 ; .Marigo et al.1 120031 ; lOhkubo et al l 
2008: (Bahena fc Klappl ,20101 ). In this paper, we focused on 
the vibrationa l stability of their evolution, and extended the 
iBaraffe et alt s analysis toward the more massive stars. That 
is, we constructed the stellar models with ^ 500Mq and car- 
ried out the linear nonadiabatic analysis against radial and 
nonradial oscillations. 
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Table 1. Property of zero-metallicity ZAMS models. 
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M cc denotes convective core mass. 
2 EQUILIBRIUM MODELS 

We adopted a Henyey-ty pe c ode developed by 
lUmeda fe Nomotol l|2002l . 120051 ) and lOhkubo et all (12009 1 
to calculate the stellar evolution. Wc constructed zero- 
metallicity models with 500, 1000 and 3000 Mq during the 
core hydrogen-burning stage. These stars were calculated 
without considering mass loss. Overshooting is not taken 
into account in this study. 

For the zero-metallicity stars, the only possible way of 
starting hydrogen burning is the pp-chain no matter how 
massive the star is. Hence, the temperature near the stel- 
lar centre is much higher than in the case of metal-riched 
stars. During the pre-main-sequence stage, gravitational 
contraction releases energy so that the energy equilibrium 
in the whole star can be kept. When the central tempera- 
ture reaches ~ 10 7 K, the pp-chain is activated and becomes 
to release enough energy to stop the gravitational contrac- 
tion for stars with < 20Mq. For the more massive stars, 
however, the gravitational cont raction cannot be stopped 
only by the pp-chain burning |Marigo et al.l l2~00ll ). When 
the central temperature reaches ~ 10 8 K, the triple alpha 
reaction produces enough 12 C to activate the CNO cycle and 
then the contraction stops. Thus, the main energy source for 
core hydrogen-burning of population III massive stars is the 
CNO cycle. Table [1] shows property of the equilibrium mod- 
els for the ZAMS stage, which we define as the time when 
the gravitational contraction stops. The central temperature 
is hardly different among the different mass models and cor- 
reponds to the occurence of the triple alpha reaction. Due 
to the large nuclear energy generation, the convective core 
is very large and occupies ~ 90 per cent of the stellar mass. 
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Figure 1. Top: variation of the growth rate of the radial fun- 
damental mode for the 500, 1000 and 3000Mq stars. The open 
triangles correspond to the ZAMS models, shown in table[T] while 
the open circles correspond to the most unstable ones, for which 
the property of oscillations is shown in table(2] Bottom: variation 
of the nuclear energy generation rate at the centre e c and the 
ratio of the central density to the average density in the whole 
star p c /(p) for 500M Q . 



part of a, ctr, represents the oscillation frequency, and the 
imaginary part of a, o~i, gives the growth rate or the damp- 
ing rate, dependi ng on its sign. We used a Henyey-type 
code developed by ISonoi fe Shibahashl (|2012h to solve this 
eigenvalue problem. Note that for stability analysis related 
to the e-mechanism, it is important to evaluate the tem- 
perature and density dependences of the nuclear reactions 
in the oscillation time sc ale. The detail is introduced by 
Sonoi fe Shibahashl (j2012T ). 

We adopted the "frozen-in convection" approximation, 
in which perturbations of the convective flux are neglected. 
The very massive stars have the large convective core and 
the convection-pulsation interaction might affect the vibra- 
tional stability. But the uncertainty in the physical process 
remains still now, and hence we neglected it for simplicity. 



3 LINEAR NONADIABATIC ANALYSIS OF 
RADIAL /NONRADIAL OSCILLATIONS 

We carried out a linear nonadiabatic analysis of radial 
(I — 0) and nonradial (I = 1, 2) oscillations. Any eigen- 
function of linear oscillations of stars is expressed in terms 
of a combination of a spatitial function and a time varying 
function. The latter is expressed by exp(iat), where a de- 
notes the eigenfrequency. The former, the spatial part, is a 
radial function for the radial oscillations, while decomposed 
into a spherical harmonic function, which is a function of the 
colatitude and the azimuthal angle, and a radial function for 
the nonradial oscillations. 

For the radial (non-radial) oscillations, equations gov- 
erning the radial functions lead to a set of fourth-order 
(sixth-order) differential equation, of which coefficients are 
complex, including terms of frequency a. Together with 
proper boundary conditions, this set of equations forms a 
complex eigenvalue problem with an eigenvalue a. The real 



4 RESULTS 

4.1 Variation of stability with stellar evolution 

As described above, the gravitational contraction continues 
while only the pp-chain works. When the CNO cycle is ac- 
tivated, the contraction stops since the energy generation of 
the CNO cycle is much higher than the pp-chain and can 
maintain the energy equilibrium without the contraction. 

Such high energy generation induces the vibrational in- 
stability due to the e-mechanism; the nuclear energy gener- 
ation increases at shrinking phase in pulsations. This makes 
expansion in the pulsations stronger and the nuclear energy 
generation drops more. This cyclical behaviour makes the 
star act as a heat engine and the pulsation amplitude will 
grow. 

Fig.[T]shows the variation of the growth rate, —ai/(2n), 
of the radial fundamental mode for 500, 1000 and 3000Mq 
stars and the profiles of the equilibrium models for 500Mq 
star with the stellar evolution. The positive value of the 
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Table 2. Property of oscillations for the models marked with the open circles in figure [T] 
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Figure 2. Workintegral W for the radial fundamental mode (top) 
and nuclear energy generation rate e (bottom) for 500Mq. The 
workintegral is normalized with the oscillation energy of the mode 
in the whole star E. 



growth rate means instability, while the negative value sta- 
bility. As the nuclear energy generation abruptly increases 
just before the ZAMS stage, the growth rate does as well 
and becomes positive. 

Fig.[2]shows the work integral for the fundamental mode 
and the nuclear energy generation rate for 500Mq. Before 
the abrupt increase in the nuclear energy generation rate, 
the e-mechanism does not work efficiently in the core. But 
after the increase, excitation by the e-mechanism causes the 
instability. 

Table[2]shows the results of the stability analysis for the 
models marked with the open circles in Fig. [T] The growth 
time-scale is defined as r g = —2n/ai, that is, the e-folding 
time of the amplitude. The positive values of the growth 
time-scale mean the instability as well as the growth rate. 
The radial fundamental mode becomes unstable since the 
amplitude is relatively large in the core and excitation by the 
e-mechanism there exceeds flux dissipation in the envelope. 
On the other hand, the other radial modes are vibrationally 
stable since the amplitude is much larger in the envelope 
than in the core (Fig. [3}. 

Although non-radial modes of I = 1 and 2 were also an- 
alyzed, no unstable modes were found. The acoustic waves of 
the non-radial p-modes can propagate in a less deep interior 
than those of the radial oscillations since they are reflected 
at the depth where the oscillation frequency is equal to the 




Figure 3. Profiles of de^/dr = a 2 ^pr 2 , which is proportional to 

the local kinetic energy density of the oscillations and normalized 
rR 

as / de^/dr dr/R = 1, of radial (I = 0) and non-radial (I = 1, 2) 
modes for the 500AfQ model marked with the open circle in Fig. 
[T] The dotted vertical line denotes the top of the convective core. 



Lamb frequency. The gravity waves of the g-modes cannot 
propagate in the central nuclear burning region due to the 
large convective core (Fig. [3)l . 

The instability of the fundamental mode continues dur- 
ing the early stage of the core hydrogen burning. The most 
unstable models, marked with the open circles in Fig. [1] are 
at age ~ 5 x 10 5 yr and the gr o wth t ime-scale is ~ 10 4 yr. As 
pointed out by Bar affe et al.l (|200ll ). the growth time-scale 
in the case of the e-mechanism is relatively long. And since 
the central temperature does not increase efficiently from 
500 to 3000Mq , the effect of the nuclear energy generation 
and the growth rate does either. 

After the activation of the CNO cycle, the stellar ra- 
dius expands. With this, the density contrast between the 
inner and outer parts of the star, or the ratio of the central 
density to the average density in the whole star, p c /(p), in- 
creases. This causes the amplitude in the envelope becomes 
much larger than in the core. Fig. [2] shows that in such a 
stage flux dissipation in the envelope exceeds the excitation 
by the e-mechanism, which works efficiently in the core. The 
instability disappears at age ~ 1 x 10 6 yr, or X c ~ 0.4. Ac- 
cording to Bar affe et all the duration of the unstable phase 
becomes longer as the stellar mass increases from 120 to 
500Mq. In this mass range, the stellar structure may be- 
come more favourable to the e-mechanism instability as the 
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stellar mass increases. In the mass range analyzed in this 
study, on the other hand, the stellar structure and hence 
the duration of the unstable phase do not change substan- 
tially with stellar mass. 

Compared with the result of iBaraffe et all the dura- 
tion of the unstable phase for 500Mq is longer, while the 
maximum value of the growth rate is lower. In their anal- 
ysis, around the ZAMS stage, the growth rate initially has 
a positive and extremely high value, and then shows the 
rapid variation with time due to the adjustment of the stel- 
lar structure when the initial CNO is produced (see their 
fig. 4). In our analysis, on the other hand, the growth rate is 
initially negative and becomes positive just after the onset 
of the CNO cycle. This implies a difference in the methods 
to calculate evolutionary models. 

4.2 Mass loss 

lAppenzellerl (|l970al lbT) performed non-linear analysis of vi- 
brational instability of very massive main-sequence stars 
against radial pulsation. It was demonstrated that recur- 
ring shock waves in the outer layers results in the forming of 
a continuously expanding shell around the star when the ve- 
locity amplitude exceeds the sound velocity. Eventually the 
outermost layers of the shell reach the escape velocity and 
mass loss occurs. In addition, due to the shock dissipation, 
the amplitude cannot increase significantly above the value 
at which the shell is form ed. 

IBaraffe et all l^OPl! ) performed linear stability analy- 
sis of the population III very massive stars and estimated 
mass loss due to the vibrational instability against the ra- 
dial fundamental mode. Although the absolute value of the 
amplitude cannot be determined in the linear analysis, they 
set the velocity amplitude to the sound velocity at the stel- 
lar surface. Besides, they assumed that the pulsation energy 
gained through the heat-engine mechanism was transformed 
into the kinetic energy of the escaping matter. 

Although non-linear analysis is actually necessary to in- 
vestigate the pulsational mass l oss, we estima ted the mass- 
loss rate in the same method as IBaraffe et al.l The mass loss 
rate was evaluated for each unstable model and integrated 
with time. In this study, the equilibrium models evolve with- 
out mass loss although their mass should be actually reduced 
by following the mass loss rate at each evolutionary stage. 

The results are shown in table [3] AAf denotes the total 
mass lost during the unstable stage. The one for 50 0Mf^ is 
roughly consistent with the result of iBaraffe et all 25Mq. 
Besides, in their result, the ratio of the mass-loss amount to 
the total stellar mass increases from 120 to 500Mq. Extrap- 
olating from their result, significant mass loss was expected 
to occur for stars with > 500Mq. Such a tendency is also 
found in our result; 4 per cent for 500Mq, while 7.7 per cent 
for 3000A/q. Such mass loss would to some extent cause 
difference from the evolution without the mass loss, but, by 
itself, might not have significant influence on the qualitative 
later scenario of the stars. 

The very massive stars have significantly large convec- 
tive cores. If the mass loss is so strong that the stellar sur- 
face would become close to the top of the convective core, the 
CNO-elements produced by the triple alpha reaction and the 
CNO cycle will reach the surface by the convective mixing. 
This could cause the radiative mass loss to become efficient. 



Table 3. Property of mass loss due to the instability of the fun- 
damental mode. 
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From left to right the table entries read: stellar mass, duration of 
vibrationally unstable phase, mass-loss rate for the model marked 
with the open circles in Fig. [T] total amount of the mass loss and 
convective core mass at the evolutionary stage when the instabil- 
ity disappears. 

For the stars with 500 — 3000Mq , the convective core occu- 
pies ~ 90 per cent of the stellar mass at the ZAMS stage (ta- 
ble [TJ , but monotonically decreases during the core hydro- 
gen burning stage. Then, the mass fraction of the convective 
core becomes less than 70 per cent at the evolutionary stage 
when the instability disappears (table [3}. Thus the amount 
of the mass loss in our result is not enough to trigger the 
above situation. For the more massive stars, however, the 
convective core and amount of mass loss are expected to be 
larger. Besides, although the instability disappears due to in- 
crease in the density contrast as discussed in section I4TT1 the 
stronger mass loss could maintain the vibrational instability 
by limiting the density contra st, like models of Wolf-Rayet 
stars shown in iMaederl Then, it is worth analyzing 

if the above situation occurs for stars with > 3000Mq . 



5 CONCLUSIONS 

We have analyzed the vibrational stability of the main- 
sequence stage of population III 500 — 3000Mq stars against 
the radial and non-radial modes. We found that only the 
radial fundamental mode becomes unstable due to the e- 
mechanism. The instability appears just after the pre-main- 
sequence contraction stops and the nuclear energy genera- 
tion becomes large enough to maintain the energy equilib- 
rium, while it disappears in the middle part of the core hy- 
drogen burning stage as the density contrast between the 
inner and outer parts of the stars becomes larger. The other 
radial and non-radial modes never become unstable since the 
amplitude in the envelope is much larger than in the core. 
Unfortunately, the result contains uncertainty due to the 
large convective core. Improvement in our knowledge about 
the convection-pulsation interaction is required to obtain a 
definite conclusion on the stability analysis. 

To measure the significance of the instability, we esti- 
mated the a mount of mass loss due to the instability by 
adopting the IBaraffe et al.l (|200ll ) method, although non- 
linear analysis was actually necessary. We found that the 
amount is less than 10 per cent of the whole stellar mass. 
Thus, we conclude that mass loss due to this instability, by 
itself, is n ot significant f or the later evolution. However, we 
note that iMarigo et alj (|2003l ) suggested that the radiative 
mass loss is supposed to siginificantly affect the evolution 
of the populat i on II I very massive stars with > 750Mq. 
lEkstrom et alj (|200fet ) showed that rotational models of pop- 
ulation III stars favour the surface enrichment of the heavy 
elements due to rotational mixing processes, which might 
trigger stronger radiative mass loss. Then, the effect of the 
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vibrational instability should be combined with the above 
effects for more exact investigation of the evolution. 
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